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Abstract. The interaction between two-level systems (TLS) and strain fields in a solid is 
contained in the diagonal matrix element of the interaction hamiltonian, S, which, in general, 
has the expression 8 = 2[y] : [S], with the tensor [7] describing the TLS "deformability" and [S] 
being the symmetric strain tensor. We construct [7] on very general grounds, by associating to 
the TLS two objects: a direction, t, and a forth rank tensor of coupling constants, [[R]]. Based 
on the method of construction and on the invariance of the expression of 5 with respect to the 
symmetry transformation of the solid, we conclude that [[R]] has the same structure as the tensor 
of stiffness constants, [[c]], from elasticity theory. In particular, if the solid is isotropic, [[R]] has 
only two independent parameters, which are the equivalent of the Lame constants. Employing 
this model we calculate the absorption and emission rates of phonons on TLSs and show that 
in isotropic solids, on average, the longitudinal phonons interact stronger with the TLSs than 
the transversal ones, as it is observed in experiments. We also show that in isotropic solids, a 
transversal wave leaves unperturbed all the TLSs with the direction contained in one of the two 
planes that are perpendicular either to the wave propagation direction or to the polarization 
direction and that a longitudinal strain applied to the solid polarises the TLS ensemble. 



In a temperature range around ten Kelvins or below, the physical properties of amorphous 
materials differ significantly from the properties of crystals and show striking universal features 
(see [l| for a collection of reviews). Most of these features can be explained by assuming that in 
the amorphous solid exists a collection of dynamic defects, which are atoms or groups of atoms, 
oscillating in two well potentials. At low temperatures the thermal activation is suppressed 
and the oscillation happens by quantum tunnelling from one potential minimum to the other, 
forming in this way what is called a two-level system (TLS). In the two-dimensional Hilbert 
space spanned by the ground states of the two wells, the effective Hamiltonian of the TLS is 





where a x and o~ z are Pauli matrices, while A and A are called the asymmetry of the potential and 
the tunnel splitting, respectively. The Hamiltonian (TjQ) may be diagonalized by an orthogonal 



transformation O, H' TLS = O t HtlsO = |cr 2 , to obtain the excitation energy of the TLS, 
e = V A 2 + A 2 -by the superscript T we denote in general the transpose of a matrix. The 
parameters A and A do not have the same values for all the TLSs, but are distributed with 
the density VP(A, A), where V is the volume of the solid. According to the standard tunneling 
model (STM), this distribution is P(A,A) = Po/A, with Pq a constant. If expressed through 
the variables e and u = A/e, the distribution function becomes P(e,u) = Pq/(u\/1 — u 2 ). 

A phonon, or any other strain in the solid body, perturbes Htls by H\ = (5/2)a z 
0; S S; 0; H: The perturbation 8 is linear in the strain field, Sy, [lj; 0] and in general may be 
written as 5 = 27^ SV,— here, as everywhere in this paper, we assume summation over repeated 
indices. The 3x3 symmetric strain tensor is defined as Sy = ^(diUj + djUi), with U{ (i = 1, 2, 3) 
being the components of the displacement field. In dyadic notations, 5 = 2 [7] : [S], where by [•] 
we shall denote matrices or second rank tensors. 

In the construction of [7] we have to be more careful. In the STM, if H\ is used to describe 
the interaction of TLS with phonons in three-dimensional (3D) bulk systems, the tensor [7] is 
replaced by a scalar, 7^ or jt, depending whether the phonon has longitudinal or transversal 
polarization, respectively, while S takes the value of the amplitude of the strain field. So one 
would write 5 = 2{^/iSi + 7tSt). 

This picture cannot be applied to the interaction of TLS with arbitrary strain fields, at 
least because simple coordinate transformations may lead to ambiguity. It is well known from 
elasticity theory that by coordinate transformations longitudinal stress can be transformed into 
shear stress and vice-versa. To avoid this ambiguity, we have to keep 7 in the form of a tensor, 
not a simple matrix, and try to find its components [6]. So, the question we ask ourselves 
is how can we extract, on very general grounds, the tensor of interaction constants from the 
physical characteristics of the TLS-note that a tensor changes its components at the coordinate 
transformations while the value of S should be the same in any coordinate system. Where can 
we begin to construct such a tensor? 

Following the line of arguments from 0] , we start by noting that since the TLS is represented 
as a tunnelling entity between two potential wells-the transition from one classical equilibrium 
position to the other taking place either by a rotation or by a translation-we may associate 
to the TLS a direction in space, t. So we have now three components, t\, t2, £3, which 
change under coordinates transformations. The simplest 3x3 symmetric tensor that we 
can construct from t is [T] = t • t T (Xy = t{tj), while a general one would have the 
components 7/^ = Ri^Tiy In abbreviated subscript notations (see for example 0] for 
details), T = (t x ,t 2 ,t 2 ,2tyt z ,2t z t x ,2t x t y ) T , and Rijki becomes Ru in a straightforward way; 
in these notations [7] becomes the vector 7 = [R] T • T. Since [S] is also transformed into 
S = (S xx , S yy , S zz ,2S yz , 2 S zx , 2S xy ) T , the component of the interaction hamiltonian, 5, is written 
simply as 5 = 2T T ■ [R] ■ S, which, we say it again, is a scalar under coordinates transformations. 

Now notice the analogy between 5 and the elastic energy density, u, that exists in a deformed 
body. In abbreviated subscript notations, u = ^S T ■ [c] • S, where [c] is the 6x6 matrix of the 
elastic stiffness constants. Under a coordinate transformation, S transforms into S' = [N] ■ S, 
where [N] is the 6x6 matrix of the transformation (see Ref. [7|, Eq. 3.34), so 

u = u' = i(S') T • [c] • S' = \S T ■ [Nf ■ [c] ■ [N] • S. (2) 

Since ([2]) should be valid for any S and any transformation, then [c] should remain unchanged- 
[c] = [N] T ■ [c] ■ [A^J-under the symmetry transformations of the crystalline lattice. From this 
argument follow all the properties of the matrix [c] which are characteristic to the symmetries 
of the crystal under consideration [3] . 

The same is true for 5. Here T transforms in the same way as S under coordinates 
transformations-T' = [N] ■ T. Therefore 5 = 5' = T T ■ [N] T ■ [R] ■ [N] • S and, like for [c], 



we have [R] = [N] T ■ [R] ■ [N] for any symmetry transformation [N]. In conclusion, [R] has the 
same structure as the tensor of elastic stiffness constants (see also [6J ) . 

If a material is isotropic, [c] has only two independent parameters-the Lame constants, A and 
/j,: cu = 2//<5/j + A for J, J < 3 and c/j = <5/j// for I or J bigger than 3. In normal subscripts 
Cijki can be written in the more compact form, Cijki = Xd~ijb~ki + ^{Sik^ji + Su^jk)- Based on the 
arguments above, the same is true for [R\; let us denote the independent parameters of [R] by £ 
and £, and, in normal subscripts, Rijki = C^ij^ki + £,(5ik$ji + &u&jk)- For a more straightforward 
reference to the STM, it might be even more convenient to denote 7 = £ + 2£ and introduce the 
reduced tensor [[r]] by r ijkl = R ijk ih = C,'$ijhl + £,'{kk&ji + kl&jk), whe I, e C + 2 C = 1 - 

Using such a form of the coupling constants tensor, we calculated in [6J the scattering rate of 
phonons on the TLSs and showed that we recover the results of the STM after we average over 
the directions of the TLS, assuming that they are isotropically oriented. The two constants of 
the STM, 7; and 74, are related to the parameters in this model by 7; = 7C; and 74 = 764, with 

Q = -Ul5 - 40£' + 32(0 2 ) and C t = ^'f . (3) 
15 15 

From Eqs. ^ we see immediately that C\ > Ct > for any real £, as it is observed 
experimentally in general 0]. 

What is also interesting to note, is that by calculating 7^ and jt from the experimental data, 
we can calculate £ and £ and determine completely the tensor of coupling constants, [R]. For 
example the two different sets of values for Po7z and Po7z in fused silica, reported by Golding et 
al (Polf = 1.4 x 1(T 5 J/m 3 and P 7 t 2 = 0.63 x 1(T 5 J/m 3 ) and Hunklinger and Arnold § 
{Poll = 2-0 x 10~ 5 J/m 3 and Po7? = 0.89 x 1CT 5 J/m 3 ), cited also by Black give the same 
solutions for £^ = 0.55 and £,' 2 = 

If we now go back to the interaction of a single TLS with a strain field and we assume that 
a transversal wave of strain field S = (0, 0, 0, S, 0, 0) T , is propagating through the solid, we find 
5 = Aj^'tytzS. This means that any TLS of t y = or t z = (i.e. any TLS which is contained in 
a plane perpendicular either to the propagation direction or to the polarization direction) will 
not be perturbed by this wave. 

Observing that for both values of £' calculated above, the corresponding values of £'(= 1 — 2£') 
are negative, we show an interesting polarization effect of the TLS ensemble, due to an 
external stress. Let us assume that we apply a longitudinal stress along the z direction, 
S = (0,0,5,0,0,0) T . This stress gives a perturbation 5 = 2^t 2 z St 2 z + 27C'5(i^ + t 2 y ) and we 
notice that the two terms in the expression of 8—5\ = 2^/t 2 St 2 and 62 = 22^/( i '(t 2 , + ty)S— have 
opposite signs. This means that, if e.g. jt 2 S > 0, the energy splitting, and therefore the 
excitation energy, of the TLSs oriented along the strain increase, while the energy splitting and 
the excitation energy of the TLSs oriented perpendicular to the strain direction decrease. In 
other words, the strain polarises the TLS ensemble. 

In conclusion, we used a model for the interaction of two-level systems (TLS) with arbitrary 
strain fields, introduced in 0], which assumes that to any TLS is associated a direction, t, and 
a tensor of coupling constants to the strain field, [[R\], and we showed on general grounds that 
[[R\] has the same structure with respect to the symmetry transformations of the solid as the 
tensor of stiffness constants, [[c]], from the elasticity theory. Some immediate consequences of 
this formalism are that in isotropic solids, on average, the longitudinal phonons interact stronger 
with the TLSs than the transversal ones, (7/ > 7t, in the language of the standard tunnelling 
model), a transversal wave does not interact with the TLSs contained in one of the two planes 
that are perpendicular either to the wave propagation direction or to the polarization direction, 
and a strain applied to the body may polarize the TLS ensemble. 
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